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ABSTRACT 

Context. Previous Spitzer and IRAS observations of the LMC suggest an excess of FIR emission with respect to the gas surface density traced by 
12 CO rotational emission lines and Hi 21 cm emission. This so-called "FIR excess" is especially noticeable near molecular clouds in the LMC, 
and has usually been interpreted as indicating the presence of a self-shielded H 2 component not traced by CO in the envelopes of molecular clouds. 
Aims. Based on Herschel HERITAGE observations taken as part of the Science Demonstration Phase, we examine the correlation between gas 
and dust surface densities at higher resolution than previously achieved. We consider three additional possible causes for the FIR excess: X factor, 
FIR dust emissivity, and gas-to-dust ratio variations between the diffuse and dense phases of the ISM. 

Methods. We examine the structure of NT80 and NT71, two molecular clouds detected in the NANTEN 12 CO survey of the LMC. Dust surface 
density maps were derived from the HERITAGE data. The gas phase is traced by MAGMA 12 CO and ATCA+Parkes Hi 21 cm observations of 
the LMC. These data provide unprecedented resolution (1') to examine the structure of molecular clouds. The dust emissivity, gas-to-dust ratio, 
and X factor required to match the dust and gas surface densities are derived, and their correlations with the dust surface density are examined. 
Results. We show that the dust surface density is spatially correlated with the atomic and molecular gas phases. The dust temperature is consistently 
lower in the dense phase of the ISM than in the diffuse phase. We confirm variations in the ratio of FIR emission to gas surface density derived 
from Hi and CO observations. There is an excess of FIR emission, spatially correlated with regions of intermediate Hi and dust surface densities 
(Ay = 1-2), and little or no CO. While there is no significant trend in the dust emissivity or gas-to-dust ratio with dust surface density, the X factor 
is enhanced at A v = 1-2. We conclude that H 2 envelopes not traced by CO and X factor variations close to the CO boundary may be more likely 
to cause these deviations between FIR emission and gas surface density than gas-to-dust ratio or emissivity variations. 
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Dust, neutral atomic hydrogen (Hi), and molecular hydrogen 
(H2) are the prime constituents of the interstellar medium in 
galaxies out of which stars form, but their amounts are often 
poorly known. In dense clouds, dust shields both H2 and its 
tracer CO from dissociation by the ambient interstellar radiation 
field (ISRF). Unlike CO, H2 is also strongly self- shielding. 
In the solar neighb orhood, H2 forms at Ay > 0.14, while CO 
requires A v > 0.8 (IWolfire et alJ2010h . Molecular clouds (MCs) 
thus consist of dense cores where CO and H2 coexist and less 
dense envelopes of H2 with little or no CO. In lower-metallicity 
environments with strong irradiation, the poorly shielded CO 
fills a much smaller fraction of the H2 volume. In those galaxies 
use of a standard conversion factor Xco to estimate H2 column 
densities from observed CO emission ca uses large amounts of 
H2 to be missed (see e.g jGlover & Mac Lowl feOlO). 

The nearest low-metallicity galaxies are the Large 
Magellanic Cloud (LMC) and the Small Magellanic Cloud 
(SMC) with [C] and [O] abundanc es 0.2 5/0.50 and 0.10/0.25 
relative to solar abundances (Pagel 2003) and distances of 50 
kpc (ISchaefedl2008b and 62 kpclszewczv k et al.ll2009h respec- 
tively. The FIR emission from dust has been used to establish 
that indeed much H2 is not traced by CO and exhibits a so-called 
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"FIR excess", implying X factors 3-6 (LMC) and 20-60 ( SMC) 
times higher tha n in the solar neighborhoo d (llsraell 1 19971; 
lLerovetal.1 120071 120091: iBernard et al.1 120081) . Their analysis 
assumes that dust grain emissivity and gas-to-dust ratios are the 
same in dense H2 clouds and more tenuous Hi regions. Our 
goal in this Letter is to explore whether these assumptions are 
justified. To this end, we ex amine the structu re of two MCs in 
the LMC, NT80 an dNT71 (Fukui et alJl2008h . Both clouds are 
relatively quiescent (Ka wamura et al.l l2009), with star formation 
rates implied by Ha and 24 /mi emission of 0.018 and 0.042 
M©/kpc 2 /yr. NT80 is located in a direction practically devoid of 
Ha emission, while NT7 1 is associated with the faint filamen- 
tary Ha nebula DEM 110 dDavies et al.ll 19761) . We examine the 
correlation between dust and gas based on HERschel Inventory 
of The Agents of Galactic Evolution (HERITAGE) data and 
MAGellanic Mopra Assessment (MAGMA, PI. T. Wong) 12 CO 
data. These data provide unprecendented resolution (15 pc) to 
observe the structure of MCs in the LMC. 



2. Observations and dust and gas surface densities 

The dust surface density (Zdust) and tempera ture Td us t in the 
LMC have been derived in Gord on et al.l (l2010h . based on Spitzer 
MIPS 160 pm observations from th e Surveying the Age nts of 
Galactic Evolution project (SAGE, iMeixner et alJ 120061) . and 
more recent SPIRE (iGriffin etluH 1201 Oh observations taken by 
Herschel as part of the HERITAGE key project during the 
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Science Demonstration Phase dMeixner et al.l 12010). The dust 
temperature was obtained by fitting a modified black body of 
emissivity law p = 1.5 to the MIPS 160 /mi, SPIRE 250 and 
350 /mi bands (the 500 /mi band being affected by an excess 
of unknown origin). The dust surface density was derived from 
the MIPS 160 /mi brightness and the dust temperature, assuming 
that the grains are silicates of density 3 g cm" 3 , size a = 0.1 /mi, 
and emissivity at 160 /mi e® 60 = 1.7xl0" 17 m 2 (absorption effi- 
ciency <2i6o = 5.47xl0" 4 ). Figures [T] and [2] show the dust surface 
density and temperature for NT80 and NT7 1 . 
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Fig. 1. Dust surface density (left), ratio of the gas surface density 
implied from dust measurements, Sjjjjf , to the gas surface den- 
sity derived from Hi and CO observations, L°^ (middle), and 
dust temperature (right) for NT80. The 3<x level in dust surface 
density is 0.04 M©/pc 2 . The purple contours show the 20, 30, 
and 40 M©/pc 2 H i surface density. The green contours show the 
10, 40, and 80 M©/pc 2 H2 surface density inferred from CO. The 
solid black lines show the MAGMA coverage. 
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Fig. 2. Same as Fig. [T] for NT71, except that the purple contours 
represent the 15, 25, and 35 M /pc 2 levels of Hi surface density. 
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Fig. 3. Gas surface densities of the atomic and molecular phases 
vs dust surface density. The dashed line represents the gas-to- 
dust ratio. 



The H i column density was taken from the Australian tele- 
scope compact array (ATCA)+Parkes map of the LMC by 
Kim et al. (2003), and converted into a surface density via X(H i) 
= 1.08xl0" 20 N(Hi), where 2(Hi) is the Hi surface density in 
M©/pc 2 , and the conversion includes the contribution of He to 
the mean molecular weight (1.36). We applied the same back- 
ground subtraction to the dust and H i surface density maps to 
set the zero level of the sky background at the end poi nts of the 
HER ITAGE scans, located outside of the LMC (Mei xner et al.l 
l2010h . The molecular gas surface density was derived from 
MAGMA CO observations via £(H£°) = 2.16xlO- 20 X co Ico, 
where ^(H^ ) is the molecular gas surface density in M /pc 2 , 
Ico is the CO integrated intensity in K km/s, and Xco is the X 
factor. We assume Xco values derived from a virial analysis of 
NT80andNT71 bv Hug hes etaPd2010h : X co = (5.1±0.1)xl0 20 
cm" 2 K" 1 km- 1 sforNT80and(4.1+0.1)xl0 2Q cm - 2 K" 1 km" 1 s 
for NT7 1 , consistent with the range of values from llsrael dl997l) 
for MCs similar to NT80 and NT71. The implied gas surface 
density is E^ s s = £(H i) + E(H£°). The sensitivities of the H i and 
MAGMA maps are 0.9 and 5.5 M©/pc 2 (0.5 K km/s). Atomic and 
molecular gas surface densities are shown in Figs. Q] and [2] 
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Fig. 4. Correlation between dust surface density and temperature 
for NT80 (black) and NT71 (red). 



3. Dust/gas correlation and FIR excess 

The first panels of Figs. [T] and [2] show that the molecular phase 
traced by CO is very well correlated with the highest dust sur- 
face density regions. The Edust > 0-08 M /pc 2 (Ay > 0.8) contour 
is indeed almost identical to the 5.5 M©/pc 2 contour of ^(H^ ) 
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NT80 



NT7' 




Fig. 5. Correlations between the dust surface density and i) the X 
factor (top), dust emissivity (middle), and gas-to-dust ratio (bot- 
tom) in NT80 (left), and NT71 (right). The dashed lines were 
obtained from the mean gas surface density, the red curves from 
the binned trends. Error bars show uncertainty in the mean. The 
shaded areas correspond to Edust < O.O5M /pc 2 , where zero 
level offsets between H i and FIR dominate. Assumptions made 
are indicated in the corresponding panels. 



(sensitivity limit). This spatial correlation is expected from the 
physics of CO formation and dissociation. The H i envelope of 
the clouds is more extended than the CO regions, but is also well 
spatially correlated with the dust surface density. Figure [3] shows 
the pixel-to-pixel correlation between the dust surface density, 
the Hi surface density, S(Hi ), the H2 surface density derived 
from CO observations, S(H^ °), and the total gas surface density, 
Xg^. £(H 1 ) dominates the gas surface density and increases lin- 
early with Sdust for Sdust < 0.1 M Q /pc 2 (Ay = 1), at which point 
the gas surface density becomes dominated by H2. The total gas 
surface density is linearly correlated with the dust surface den- 
sity over the entire range of dust surface densities. The slope 
of the correlation gives the gas-to-dust ratio, the value of which 
is GDR = 351+5 for NT80, and GDR = 234+4 for NT71. The 
intercept of the total gas surface density with the zero dust sur- 
face density is at £^ s s = 1.4+0.56 M /pc 2 for NT80 and 3.8+0.5 
M /pc 2 for NT7 1 , indicating that there is an offset between the 
zero levels of the dust and H 1 surface density maps. As a result, 
we do not trust ratios of dust and gas surface densities at low 
surface densities (Edust < 0.05 M /pc 2 ). 

We examine the correlation between two different estimates 
of the gas surface density: from dust measurements and a con- 
stant gas-to-dust ratio, Ejjjjf = GDRx£ du st, and from co and 
Hi observations, E°!^. The middle panels of Figs. [T] and [2] show 

the ratio EjjSf/E^. On average, the ratio 1^/1^ is one, with 
some deviations that appear spatially correlated with the differ- 
ent phases of the ISM. In particular, EJ^VSw is highest (>1.5) 



in regions with intermediate dust (Sdust ~ 0.1-0.2 M©/pc 2 ) and 
Hi (2(Hi) = 20-30 M©/pc 2 ) surface densities, and little or no 
CO (E(H£°) < 10 M Q /pc 2 ). It is close to one (> 0.7 and < 1.3) 
at high dust surface densities (Sdust > 0.2 M©/pc 2 or Ay > 2), in- 
side the CO boundary (E(H£°) > 10 M Q /pc 2 ). It is low (< 0.5) in 
diffuse regions, outside of the H 1 and CO contours in Figs.[T]and 
[2] A low Egaf/^gas ra ti° at l° w dust surface densities is uncertain 
as it is likely dominated by small offsets between the H 1 and the 
dust surface density zero levels. On the other hand, the excess 
of FIR emission (i.e., of dust surface density) in regions with in- 
termediate dust surface density and little or no CO supports the 
presence of H2 envelopes not traced by CO, and is consistent 
with previou s conclusions drawn from the compa rison between 
dust and gas dlsraelll 19971: ILeroy et al.ll2007il2009b . 

Last, the right panels of Figs. [T] and [2] show that the dust tem- 
perature appears to be spatially anti-correlated with the dust sur- 
face density, the high dust surface density regions being colder 
than the low dust surface density regions by a few K. This effect 
is further seen in Fig.|U which shows the pixel-to-pixel correla- 
tion between Edust an d Tdust- This anti-correlation suggests that 
the regions of MCs that are well shielded from the ambient radi- 
ation are colder than the envelopes of the clouds, more exposed 
to the ISRF. This effect has not been observed at 4' resolution in 
the dust properties derived from I RAC, MIPS, a nd IRIS obser- 
vations of NT80 and NT71 dParadis et alJl2010h . but is clearly 
seen at V resolution in our Herschel data. The dust temperature 
in NT7 1 is higher than in NT80, which may result from heating 
by star forming regions embedded in NT7 1 . 



4. Possible causes of the variations of Sdust/^! 

4.1. X factor variations 



obs 
gas 



The molecular gas surface density derived from CO observa- 
tions was computed with a constant X factor. In reality, the 
CO/H2 abundance is sensitive to photo-disso ciation at Ay < 2-3 
(Rubi o et al.lll993UGlover & Mac Lowll201Qh . As a result, the X 
factor is expected to decrease (the CO/H2 abundance to increase) 
with dust surface density in the transition region between the H2 
envelopes and the CO cores of MCs. While H2 gas not traced by 
CO in the envelopes of MCs might account for the excess of FIR 
emission with respect to the gas surface density outside the CO 
boundary, unaccounted for X factor variations may also cause 
deviations in the dust/gas correlation inside the CO boundary. 

Within the CO boundary (where I C o is above the MAGMA 
sensitivity), we derive the X factor required to match the gas 
surface density inferred from dust and a constant GDR with the 
surface density implied by CO and H 1 observations: 



Xgg = (GDR x Ed U st - 1{HI)) I (2.16 x 10- 20 /co) 



(1) 



The top row of Fig. [5] shows the correlation between X£ ^ and 
Xdrnt- The red curve indicates the binned trend (0.02 M©/pc 2 
bins). Since i) we assume a constant X factor, Xco = 5.1xl0 20 
for NT80 and X co = 4.1xl0 20 for NT71, to derive GDR, and 
ii) we assume GDR to derive X°^ at each pixel, it follows 
that the average of X°^ and the assumed Xco must be and 
are equal within the error bars (<X°^> = (5.7±0.5)xl0 20 and 
(3.6±0.5)xl0 20 for NT80 and NT71). Thus, Fig. merely in- 
vestigates whether systematic variations in X°^ with Edust can 
explain the variations in Figs. [T] [2] and the scatter in Fig. [3] 

Xco is higher in the range Edust = 0.1-0.2 M©/pc 2 (Ay = 1- 
2) by a factor of up to 8 compared to the densest regions, well 
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inside the CO boundary (Sdust > 0.2 M©/pc 2 or Ay > 2). This 
enhancement however only appears marginally significant in the 
binned trends. Nonetheless, X factor variations may very well 
contribute to the observed variations in the FIR emission/gas sur- 
face density ratio inside the CO boundary. In fact, this increase 
in Xco at intermediate surface densities is likely coincident with 
the transition regions between dissociated and shielded CO, and 
supports the presence of H2 envelopes not traced by CO. The 
decrease in Xco at low (< 0.05 M©/pc 2 ) dust surface densities is 
likely due to small offsets between the H 1 and dust surface den- 
sity zero levels — the H 1 level being slightly higher, as shown 
by the negative values of Xco- Besides being difficult to explain 
physically, we do not trust its significance. 

4.2. Dust emissivity variations 

The dust surface density was derived assuming that the emis- 
sivity of dust does not depend on environment. An emissivity 
increase in the FIR of a factor 3 to 4 between the the diffuse and 
dense phases has however been invoked to explain the cold tem- 
peratures a nd the 60 11m emission deficit observed in the molec- 
ular phase ( Ste pnik et al]|2003h . and is expecte d from grain co- 
agulation in the dense phase of the ISM (JParadis et al.ll2009l) . In 
the Milky Way, this arg ument is supported by recent FIR and 
sub-mm observations by iParadis et al.l (120091) . 

The dust emissivity per unit mass, e^, was derived from 
matching the 160 jim emission to the surface density implied by 
CO and H 1 observations for a constant gas-to-dust ratio: 



Il60 - £°60^gasBl60(Tdust)/G/}R 



(2) 



where, /160 is the brightness observed at 160 yum, Bi6o(Tdust)> is 
the Planck function at the dust temperature Td us t and at 160 /mi, 
and GDR is a constant gas-to-dust ratio. 

The second row of Fig. [5] shows the pixel-to-pixel correlation 
as well as the binned relation between ^q/ 6 ?60 anc * ^dust, wnere 
6^ 60 is the constant emissivity assumed to derive the dust surface 
density dGordon et al.ll2010h . For both NT80 and NT71, ej*® is 
constant with Edust within the scatter. Again, we do not take the 
lowest, uncertain Sdmt points into account. While it is possible 
that trends be hidden in the scatter, our data do not seem to sup- 
port emissivity variations as a major contributor to the variations 
in the FIR emission/gas surface density correlation. Further in- 
vestigation with the full extent of the HERITAGE survey will be 
necessary to draw firmer conclusions. 

4.3. Gas-to-dust ratio variations 



Our analysis in Sects. I47T1 and 14. 21 was based on the assumption 
of a constant gas-to-dust ratio. It is possible, however, that Xco 
and the FIR dust emissivity are approximately uniform, and 
that GDR varies. In this case, the middle panel of Figs. [T] and [2] 
represents the variations in gas-to-dust ratio implied by gas and 
dust observations. Gas-to-dust ratio variations could be caused 
by dust destruction (or change of size) in shocks and intense 
ISRFs in the LMC, or by grain growth in molecular cores. 

The gas-to-dust ratio implied by dust and gas observations 
was obtained via GDR obs = E^/Edust. The plausibility of 
gas-to-dust ratio variations as a cause for deviations in the FIR 
emission/gas correlation was further tested by examining the 
correlation between GDR obs and 2dust, shown in the bottom 
row of Fig. \5\ The dashed line indicates the constructed trend 
obtained for a constant, mean gas surface density. If the lowest, 
uncertain points in Edust are omitted, the gas-to-dust ratio 



appears rather constant with Edust* within the scatter. Although 
a more complete investigation is needed to draw strong con- 
clusions, gas-to-dust ratio variations between the diffuse and 
dense phases of the ISM do not appear to contribute much to 
deviations in the FIR emission/gas correlation. 



5. Conclusion 

We have examined the correlation between dust, atomic, and 
molecular gas using HERITAGE, ATCA Hi 21 cm, and 
MAGMA CO observations of two LMC molecular clouds. The 
dust temperature appears consistently lower in the dense phase 
than in diffuse regions. The dust surface density is spatially cor- 
related with the atomic and molecular phases, making HerscheVs 
angular resolution and complete coverage of the IR SED a pow- 
erful tracer of molecular gas. We have however observed an ex- 
cess of FIR emission with respect to the gas surface density im- 
plied by CO and H 1 observations, which occurs at intermediate 
dust surface densities (0.1-0.2 M /pc 2 ), outside and close to the 
CO boundary. This likely indicates that molecular clouds are sur- 
rounded by envelopes of H2 not traced by CO. The presence of 
unaccounted for H2 envelopes is further supported by an increase 
in the X factor at intermediate dust surface densities, correspond- 
ing to the transition region between dissociated and shielded CO. 
We reviewed two alternative explanations for the FIR excess: 
variations in dust emissivity and the gas-to-dust ratio between 
the diffuse and dense phases of the ISM. We derived the dust 
emissivity and gas-to-dust ratio required to match the observa- 
tions, and examined their correlations with the dust surface den- 
sity in order to evaluate the plausibility of each hypothesis. We 
found that the dust emissivity and gas-to-dust ratio in NT7 1 and 
NT80 are constant with Edust within the scatter, and conclude that 
dust emissivity and gas-to-dust ratio variations are therefore un- 
likely to be responsible for the FIR excess observed near these 
clouds. Variations in emissivity and gas-to-dust ratio between 
the dense and diffuse ISM phases cannot be definitively ruled 
out however, due to uncertainties at low dust surface density that 
are caused by offsets in the zero levels of the H 1 and dust maps. 
In the immediate future, we will conduct a full investigation of 
all these effects using detailed modeling in combination with the 
completed HERITAGE survey of both Magellanic Clouds. 
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